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Abstract 

Understanding possible electromagnetic signatures of the merging and collapsing compact object 
is important for identifying possible sources of LIGO signal. Electromagnetic emission can be 
produced as a precursor to the merger, as a prompt emission during the collapse of a neutron star 
and at the spin-down stage of the resulting Kerr-Newman black hole. For the NS-NS mergers, 
the precursor power scales as L ^ B^gGM^sR^g/{R'^^^c), while for the NS-BH mergers, it is 
(GM/(c^i?Ns))^ times smaller. 

We demonstrate that the time evolution of the axisymmetric force-free magnetic fields can be 
expressed in terms of the hyperbolic Grad-Shafranov equation and formulate the generalization 
of the Ferraro's law of iso-rotation to time-dependent angular velocity. We find exact non-linear 
time-dependent Michel-type (split-monopole) structure of magnetospheres driven by spinning and 
collapsing neutron star in Schwarzschild geometry. 

Based on this solution, we argue that the collapse of a NS into the BH happens smoothly, without 
natural formation of current sheets or other dissipative structures on the open field lines and, thus, 
does not allow the magnetic field to become disconnected from the star and escape to infinity. 
Thus, as long as an isolated Kerr black hole can produce plasma and currents, it does not lose its 
open magnetic field lines, its magnetospheric structure evolved towards a split monopole and the 
black hole spins down clcctromagnetically (the closed field lines get absorbed by the hole). The 
"no hair theorem", which assumes that the outside medium is a vacuum, is not applicable in this 
case: highly conducting plasma introduces a topological constraint forbidding the disconnection of 
the magnetic field lines from the black hole. Eventually, a single random large scale spontaneous 
reconnection event will lead to magnetic field release, shutting down the electromagnetic black 
hole engine forever. Overall, the electromagnetic power in all the above cases is expected to be 
relatively small. 

We also discuss the nature of short Gamma Ray Bursts and suggest if the magnetic field is 
amplified to ~ 10^^ G during the merger or the core collapse, the similarity of the early afterglows 
properties of long and short GRBs can be related to the fact that in both cases a spinning black 

hole can retains magnetic field for sufficiently long time to extract a large fraction of its rotation 
energy and produce high energy emission via the internal dissipation in the wind. 



2 



I. INTRODUCTION 



Estimating possible electromagnetic signature of merging and collapsing neutron stars 
is most desirable for the gravitation waves searchers by LIGO and for identifying possible 
progenitors of short Gamma Ray Bursts. Collapse of a neutron star into black hole may 
proceed either through the accretion induced collapse (AIC) or during binary neutron star 
mergers. We expect at late stages both processes proceed along a somewhat similar path: in 
case of the merger, the two collapsing neutron stars form a transient supermassive neutron 
star which then collapses into the black hole. Both an accreting neutron star (in case of 
an AIC) and the transient supermassive neutron star are expected to be magnetized. In 
addition, in case of merging neutron stars the strong shearing of the matter may increase 
magnetic field well above the initial values. 

In case merger of compact stars the electromagnetic power can be generated as a precursor 
to the merger due to either effective friction of the neutron star magnetospheres, or due to 
purely general relativistic effect, see ^|TT] Later, and in the case of the AIC, several types of 
electromagnetic emission can be foreseen. First, the electromagnetic power in vacuum may 
be generated directly, due to the changing magnetic moment of the collapsing star ^ [2] . 
Even if the outside medium is highly conducting, electromagnetic may be generated via 
effective (resistive) disconnection of the external magnetic fields, provided that the collapse 
naturally leads to formation of narrow dissipative current structure. Second, a pulsar-like 
electromagnetic power generated by the rotation of the neutron stars and extracted via the 
magnetic field. As we argue below, as long as the black hole can produce plasma via vacuum 
breakdown, it can self-generate electric currents, retain the magnetic fields and spin-down 



electromagnetically for time periods much longer than the collapse time, see 



Conventionally, in estimating the possible electromagnetic signatures it was first assumed 
that a fraction R-^s/Rg of the initial external magnetic energy (also built-up by the collapse 
and compression of the magnetic field) is radiated away on time scale of the order of the 
collapse time [3] . Ref . [^l ISj considered radiation from accelerated changes in the magnetic 
moment during collapse, producing energy E ~ i?Q(i?G'i?Ns)^''^ (somewhat smaller than the 
energy of the magnetic field before the collapse). Along the similar lines, Ref. [2] employed 
GRMHD simulations and followed a collapse of a non-rotating neutron star into the black 
hole. 
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In our view the main limitation of these models is that the external medium was treated 
as a vacuum. Electrodynamically, vacuum is a highly resistive mediums, with the impedance 
of the order of 47r/c = 477 Q. As a result, nothing prevents magnetic fields from becoming 
disconnected from the star and escaping to infinity. We expect that the magnetic field dy- 
namics would be drastically different if the external magnetosphere were treated as a highly 
conducing medium. This is a common consequence in relativistic astrophysical sources, 
since ample supply of plasma is available through vacuum breakdown. For example, inves- 
tigating the dynamics of the magnetic field in the simulations in Ref. [2] [see also E] shows 
that during the collapse the magnetic field becomes effectively disconnected from the star, 
at distances somewhat larger than the Schwarzschild radius. If the outside medium were 
treated as highly conducting plasma, such processes would be prohibited. The importance 
of resistive effects in the magnetosphere was stressed early on in the original paper by [1], 
who point out that "for spherically-symmetric collapse there is no energy released to the 
outside at all." 



Magnetic field may still escape to infinity if the collapse naturally creates conditions 
favorable for reconnection, e.g., by forming narrow current sheets or leading to the overall 
breakdown of fluid approximation by creating regions where electric field exceed magnetic 
field (the latter regions are naturally created both in pulsar magnetospheres [7] and near 
black holes moving through the external magnetic field [8]). In this paper we address a 
question " does the collapse of rotating magnetized neutron star naturally creates condition 
for efficient reconnection of magnetic field lines well before the foot points cross the horizon?" 
We argue that this does not happen. 



The plan of the paper is the following. In ^|TT] we discuss possible types of precursor 



emission in NS-NS, NS-BH and BH-BH mergers. In ^11 we make estimates of the classical 
(non-GR-modified) pulsar-like prompt electromagnetic power during collapse. In the main 
§IV| we find exact solutions for the structure of collapsing magnetospheres. Based on this 
solution we argue that as long as the resulting black hole can produce plasma and currents 
by vacuum breakdown, it may produce electromagnetic much longer that the collapse time. 



II. PRECURSOR EMISSION IN MERGERS 



For merging compact objects (NS-NS, BH-NS, BH-BH) a number of mechanism can 
generate precursor or afterglow emission. In case of merging neutron stars, one expects an 
electromagnetic precursor due to effective 'friction' of the neutron stars' magnetospheres 
P4TT]. Qualitatively, a neutron star moving through magnetic field generates an inductive 
potential drop, which induced real charges on the surface, which in turn produce a component 
of the electric field along the magnetic field and electric currents. The estimate of the 
corresponding power is 

Lu^B'^,Rl,(3\ = B'^,GM^, (1) 

where -Bns is the surface magnetic field of a neutron star, i?NS is its initial radius and M 
is its mass, P = v/c is the dimensionless velocity of a neutron star. The last equality in 
Eq. ([T]) assumes a Keplerian orbit with radius -Rorb- The estimate ([T]) can be derived by 
calculating the potential drop across the neutron star, A$ ^ /3Sns-Rns and assuming the 
resistance of the resulting electric circuit to be close to the vacuum inductance ~ 47r/c. 
Just before the contact, the unipolar power ([!]) is 

L^,^a. = 6xlO^Vgs-^(J%)' (2) 

for Mns = IAMq and i?Ns = 10 km. 

The total electromagnetic energy produced by the unipolar induction mechanism can be 
found by integrating power ([T]) with the radius evolving due to radiation of gravitational 
waves, R = -Rlc (1 — {GM)H/{c^Rf^Q))^^^ and magnetic field scaling as B = B^s{.Rns/RY 
(the model becomes applicable when the magnetospheres of the neutron stars touch at the 
light cylinder distance -Rlc! at earlier time the interaction is through winds and scales as a 
sum of the spin-down powers of the neutron stars), see [TT], 

- IZ ^"-^ - (tf )^ ^ ^ ^ 

In addition, there is a purely general relativistic effect, when the motion of the compact 
object across magnetic field in vacuum generates parallel electric field, which in turn leads 
to generation of plasma and the production of electromagnetic outflows with power [8] 



{GMfBlP^ {GMfBl 



[see also[T2l [13]. This type of interaction is important for BH-NS and BH-BH mergers, in 
which case there are no real induced charges to produce the parallel electric field, the parallel 
electric field is a pure vacuum effect, resulting from the curvature of the space-time. This 
power is smaller than for NS-NS coalescence by a factor {Rq/R-^s)'^, where Rg = 2GM/(? 
is the Schwarzschild radius. 

Qualitatively, the power Q can be estimated from the potential drop across the 
Schwarzschild horizon. There is an important difference between NS-NS and BH-NS elec- 
tromagnetic interaction, though: in case of the NS-NS system, the parallel electric field is 
produced by real surface charges [H], while in case of the black holes the parallel electric 
field is a pure vacuum effect, resulting from the curvature of the space-time [8]. 

For NS-BH system just before the contact, the general relativistic unipolar power LufiR 

is 

i:^,GiJ = 3xlO^Vgs-^(J^)' (5) 

The total emitted energy is 

(Relations ([5]j6]) assume equal masses of the merging objects; it is straightforward to gener- 
alize them to unequal masses.) Thus, the total energy dissipated via the general-relativistic 
unipolar induction mechanism is of the order of the magnetic energy of the neutron star. 
Note, that the energy is taken from the linear motion of the neutron stars, and not from the 
energy of the magnetic field. 

In addition, a more involved electromagnetic signatures are expected due to the pertur- 
bations that the merging black holes induce in the possible surrounding gas [T5H20] . 

III. PULSAR-LIKE PROMPT ELECTROMAGNETIC POWER DURING COL- 
LAPSE 

In this section we discuss pulsar-like electromagnetic power during the prompt stage of 
the neutron star collapse, treating the collapse approximately, in a classical regime up to the 
Schwarzschild radius -Rg- As the neutron star collapses, it spins up, VL oc -R~^, magnetic field 
increases due to flux conservation, B cx R~^, while the radius decreases. Let us first discuss 
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how electromagnetic power evolves during the prompt stage of the collapse, neglecting, for 
the time being, the effects of General Relativity. 

If the dipole spin-down formula remains valid, the dipolar electromagnetic power increases 
according to 

L,^B^X^' = L^s{^)' (7) 

where Lns ~ -Sns-^ns^ns/"^^ is the standard pulsar dipolar spin-down. (Note that the 
magnetic moment oc BR^ oc R decreases during the collapse.) In case of a free-fall of 
the neutron star surface, Rg = -Rns(1 — t/tc)'^^^, where tc = {2/3)R^^ / VGM, resulting in 
luminosity evolution 

Limiting the collapse to the fall time down to the Schwarzschild radius, tj = {2y/2/3){R^^ — 
R^q"^) / {c\/ Rg) 1 the total released energy is relatively small. 

i^iot — ^7^5 -t^NS ~ -t^NS-n^Ns/C \^) 

G 

The pulsar-like luminosity of the collapsing neutron star may be a bit larger than given 
by Eq. ([T]). Under the ideal MHD condition, the magnetic field is frozen into plasma. Thus, 
for field lines penetrating the star, the angular velocity of the field lines is locally equal to 
the angular velocity of the foot-point. The collapse is expected to produce strong shearing 
of the magnetic field lines' foot-points. As a result, large scale currents will be launched into 
the magnetosphere, increasing the spin-down power. Increased currents will tend to inflate 
the magnetosphere, resulting in an increased magnetic flux through the light cylinder and 
higher spin-down luminosity [21]. As the upper limit, we can use the spin-down power of 
the split monopole solution, 

_2 {BsRlfnl _( c Y{Rns\\ _{ c Y L^s f,^. 
^ 3 c ^ [r^s^^sJ \ Rs ) [r^s^^sJ il-t/Q'/' ^ ^ 

Larger currents in the magnetosphere lead to the increase of power, but for the collapsing 
neutron star the power increases slower with the decreasing radius. As a result, the total 
energy released during the collapse time tc remains fairly small 



^tot 



(11) 



--Rns^ns/ ^0 

Since the collapse time is short and R^s not much larger than Rg, the increase in power is 
mild at best, while the total released energy is small. 



IV. MAGNETOSPHERES OF COLLAPSING NEUTRON STARS 



A. Direct emission of electromagnetic waves during collapse 

As a neutron star experiences a collapse, the frozen-in magnetic field evolves with time, 
generating electric field and a possible electromagnetic signal. Historically, the first treat- 
ment of the electromagnetic fields of a collapsing neutron stars was done in the quasi-static 
approach [1] , in which case the electric field follows from the slow evolution of the magnetic 
field. The quasi-static approach was later demonstrated to give the incorrect asymptotic 
decay of the fields with time ^22j. As the neutron star contracts, the magnetic moment 
decreases oc Rg. The scaling of the decay of the fields on the BH calculated in Ref. |22] was 
confirmed by [2] , who performed numerical simulations of the neutron star collapse into BH 
and saw a predicted power-law decay of the electromagnetic fields. 

Most of the power in the calculations done in Ref. [2] was emitted at times of the order 
of the collapse time, well before the predicted asymptotic limit. Overall, the simulations are 
dominated by heavy resistivity effects intrinsic to the vacuum approximation: the discon- 
nection of the magnetic field field lines from the star typically (except in the Kerr-Schild 
coordinates) occurs at the time when the strong compression of the magnetic field against 
the horizon and the corresponding effects of the numerical resistivity becomes important. 

The assumption of a highly conducting exterior changes the overall dynamics of the 
electromagnetic fields. As we argue below, the high conductivity of the external plasma 
would prevent the formation of disconnected magnetic surfaces, formally prohibiting the 
processes described in Ref. |2j. 

B. Force-free approximation in General Relativitiy 

There is a broad range of astrophysical problems where the magnetic fields play a dom- 
inant role, controlling the dynamics of the plasma [23j. The prime examples are pulsar 
and black hole magnetospheres; Gamma-Ray bursts, AGN jet may also be magnetically 
dominated at some stage [e.g., El]. If the magnetic field energy density dominates over 
the plasma energy density, the fiuid velocity, enthalpy density and a pressure become small 
perturbations to the magnetic forces. The dynamics then can be described in a force-free 
approximation [25j. In the non-relativistic plasma the notion of force-free fields is often 
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related to the stationary configuration attained asymptotically by the system (subject to 
some boundary conditions and some constraints, e.g., conservation of helicity). This equi- 
librium is attained on time scales of the order of the Alfven crossing times. In strongly 
magnetized relativistic plasma the Alfven speed may become of the order of the speed of 
light c, so that the crossing times becomes of the order of the light travel time. But if plasma 
is moving relativistically its state is changing on the same time scale. This leads to a notion 
of dynamical force-free fields. 

MHD formulation assumes (explicitly) that the second Poincare electro-magnetic invari- 
ant E-B = and (implicitly) that the first electro- magnetic invariant is positive B^ — E^ > 0. 
This means that the electro-magnetic stress energy tensor can be diagonalized and, equiv- 
alently, that there is a reference frame where the electric field is equal to zero, the plasma 
rest frame. This assumption is important since we are interested in the limit when matter 
contribution to the stress energy tensor goes to zero; the possibility of diagonalization of 
the electro-magnetic stress energy tensor distinguishes the force-free plasma and vacuum 
electro-magnetic fields, where such diagonalization is generally not possible. 

The equations of the force-free electrodynamic can be derived from Maxwell equations 
and a constraint E-B = 0. This can be done in a general tensorial notations from the 
general relativistic MHD formulation in the limit of negligible inertia [26]. This offers an 
advantage that the system of equations may be set in the form of conservation laws [27J. A 
more practically appealing formulation involves the 3-1-1 splitting of the equations of general 
relativity [281 129] . The Maxwell equations in the Kerr metric then take the form 

V ■ E = Arrp 

V ■ 5 = 

V X (aB) = A7raj + DtE 

V X (aE) = -DtB (12) 

where Dt = dt — is the total time derivative, including Lie derivative along the velocity of 
the zero angular momentum observers (ZAMOs), V is a covariant derivative with the radial 
vector = adr and a = y^l — 2M/r. Taking the total time derivative of the constraint 
E ■ B = and eliminating DfE and DtB using Maxwell equations, one arrives at the 
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corresponding Ohm's law in Kerr metric [8], generalizing the result of [25]: 



^ i^B X (aB) - E -V X {aE)j B + a{V ■ E)E X B 

^ = ^^^^ 

Note that this expression does not contain the shift function /3. 

The generic limitation of the force-free formulation of MHD is that the evolution of 
the electromagnetic field leads, under certain conditions, to the formation of regions with 
E > B [e.g., [7], since there is no mathematical limitation on B^ — E"^ changing a sign 
under a strict force-free conditions. In practice, the particles in these regions are subject 
to rapid acceleration through E x B drift, following by a formation of pair plasma via 
various radiative effects and reduction of the electric field. Thus, regions with E > B are 
necessarily resistive. This breaks the ideal assumption and leads to the slippage of magnetic 
field lines with respect to plasma. In addition, evolution of the magnetized plasma often 
leads to formation of resistive current sheets, with the similar effect on magnetic field. If such 
processes were to happen in the magnetospheres of the collapsing neutron star, this might 
potentially lead to disconnection of the magnetic field lines form the star and a magnetic 
field-powered signal. Below we argue that in case of collapsing neutron stars this does not 
happen. 



V. THE RESTRICTED WAVE GRAD-SHAFRANOV EQUATIONS 

Let us derive a dynamic equation that describes the temporal evolution of the force-free 
fields in special relativity under the assumption that the fields remain axially-symmetric. 
Previously the equations governing general time-dependent force-free motion has been writ- 
ten by [301 E]. 



In relativistic plasma the force-free condition is given by the Ohm's law (13), where in 
this section we set a = 1. Generally, any function can be represented as a sum of a gradient 
and a curl of a vector function. Under the assumption of axial symmetry and zero divergence 
for magnetic field, we can express electric and magnetic fields as 

„ VPxe^ 21 ^ 

B = ^ 

r sm 6 r sm 6 

E = -V<1> + — ^ — 14 

r sm t7 r sm tl 
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where P is the magnetic flux function P = A^w, w = r sin 9, A^, is the electric potential and 
K and L are some arbitrary functions to be determined, / is the poloidal current through a 
flux cross-section divided by 27i. The Maxwell equation = —V x E gives 



L = dtP/2 

dtI=Ud',K+\smedg^ 
2 \ smtJ 



A* = sin^ 



"2 sin^ 9 



-A*K 
2 



(15) 



(16) 



The ideal condition E ■ B = implies 



2IdtP = - (VK + r sin x e^) ■ VP 



(17) 



Equations ( 15 )-( 17) highlight two separate types of non-stationarity: (i) due to the variations 



of the current I{t) for a given shape of the flux function (Eq. (16)); (ii) due to the variations 



of the shape of the flux function for a given current / (Eq. (17)). 



A. Constant shape of flux functions, dtP = 0, variable current 



Let us flrst consider the case when dtP = 0. Then Eq. ( 17) implies that WKo+r sin 0V$ x 
e,^ is orthogonal to VP (and is thus along the poloidal magnetic fleld). Above, Kq denotes 
a particular case when the P is constant in time. Thus 

Vi^'o = -r sin 9V^ x + r sin 9nVP x (18) 

E = -fiVP = -v^e^ X B 

T sin f) ttP' 

dtl = -^VP xVQ-e^ = ^{B -VQ) (19) 

where fl is an arbitrary function, which can be identified with the angular velocity of the 
rotation. 

The (f) component of the induction equation then becomes the time-dependent Grad- 
Shafranov equations for the restricted case when the shape of the flux surfaces remain 
constant, but the angular velocity Q and, thus, the poloidal current are time and space- 
dependent: 

.'v(i^vp) + ^^<^^ + ..^n(v/>.vn) = o (20) 
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This is a Grad-Shafranov for axisymmetric force-free structures that rotate with arbitrarily 
varying angular velocity, but keep the shape of the flux functions constant. 
The poloidal components of the induction equation give 



2VP X V/ ■ e 







ro(VP)2 



{vpy 



(21) 



Note that Eqns (19) and (21) involve only poloidal magnetic field which under assumption 
dfP = remain constant in time. 



Eqns (19) and (21) can be combined to determine the evolution of f2: 

,2 



zu 



B-V{B-Vn) 



(22) 



where Bp is the poloidal magnetic field. Eq. (22) is the generalization of the Ferraro's law 
of iso-rotation to time-dependent angular velocity. 



Eqns (19, 20, 21) constitute a closed system of equations for variables P, /, VL under the 
assumptions of time- dependent / and Vt and stationary P. Generally, it is not guaranteed 
that there is a physically meaningful solution of this system: recall that this system de- 
scribes a restricted motion of force-free plasma, when the shape of the flux function remains 
constant. 



B. Variable shape of flux functions 



By virtue of (17) and (19) variable shapes of the flux functions can be described by the 
addition to V-R'o of a term proportional to VP, K = Kq + F{P). 

Let us first consider K = F{P) separately, neglecting the cross-terms in electric field. 
The E ■ B = gives 

VP ■ VP = 2IdtP (23) 

or, since F = P(P), 

F'iWPf = 2IdtP (24) 
The Maxwell equation c?fB = —V x E gives 

dj = ^A*F = ^ (P'A*P + (VP)^P") (25) 

The (J) component of the induction equation then gives the Grad-Shafranov eq. 



(26) 
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This is a wave (hyperbolic) Grad-Shafranov for non-rotating axisymmetric force-free struc- 



tures that evolve with time. The current / here is determined from Eqns. (24)-(25). 

The wave Grad-Shafranov equation can be written in a general case, when both current 
and the flux function evolve with time (Appendix |A]), but it's overly complicated form makes 
it not useful for practical purposes. 



C. Time-dependent Michel's split-monopole solution in flat space 

Both in the case of accretion induced collapse and for NS-NS mergers, right before the 
final plunge the NS is expected to rotate with a spin close to break-up limit of ~ 1 msec. 
As a result, the light cylinder is located close to the NS surface. The theory of pulsar mag- 
netospheres predicts that outside the light cylinder the magnetic field structure resembles 
the split monopole structure [22] • This is confirmed by direct numerical simulations |33j . 

In ^|V] we derived hyperbolic wave Grad-Shafranov equation, describing time-dependent 
force-free electromagnetic fields. It may be verified directly, that the Michel's monopole so- 
lution for rotating force-free magneto sphere is valid for time- dependent angular velocity 



Vt, surface magnetic field and neutron star radius Rs- For monopole field, Eq. (22) gives 
a radially propagating fast wave 



d^Q = din 



VL = VL{r± t) (27) 
The flux conservations requires BgR^^ =const= -Bns-Rns- Then the Grad-Shafranov equation 



(20) has a slit-mo nopole-type solution for electromagnetic fields of the collapsing neutron 



star: 



Br = {^~~^ -^S' ^<i> ~ ^~^r ~ ^'f 

> = -2 f-^j cos^^]5, 



p= {1- COS 9)B,Rl 

$ = -pn 

PiP ~2BsRl)VL 1 
/ = - ^ ^^J, = -^BsRln sin' 9 (28) 

where P is the flux function, and $ is the electric potential and Q = Q{r — t). It may be 



verified directly that Eq. ( 19 ) is satisfied. 
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Thus, we found exact solutions for time-dependent non-linear relativistic force-free config- 
urations. Though the configuration is non- stationary (there is a time-dependent propagating 
wave), the form of the flux surfaces remains constant. 



VI. ELECTRODYNAMICS OF NEUTRON STAR COLLAPSE 
A. Force-free collapse in Schwarzschild metric 

Next we apply the solutions obtained in the previous section to the electrodynamics of 
neutron star collapse taking into account general relativistic effects. The split monopole 
solution may be a good approximation for several reasons. First, the collapse is likely to 
induce strong shear of the surface foot-points. As a result, strong electric current will be 
launched in the magnetosphere strongly distorting it. Highly twisted magnetic field lines will 
tent to open up to infinity, so that the magnetosphere will resemble a monopolar solution at 
each moment corresponding to the changing angular velocity of the surface foot-points. For 
a general case of strongly sheared foot-points, a time-depended angular velocity will break 
a force-balance. Still, we expect that the overall dynamical behavior will be similar to the 
time- dependent Michel's solution. 

Second, as we argue below, the open field lines cannot slip off the horizon, while the closed 
field lines will quickly be absorbed by the black hole. Thus, the magnetosphere of the black 
hole will naturally evolve towards the split monopole solution. Fig. [TJ Finally, in a more 
restricted sense, the fully analytically solvable dynamics of the monopolar magnetosphere 
collapse can be used to estimate the physical effects occurring on the open field lines. 

The stationary Michel's solution has been generalized to Schwarzschild metric [34j (BZ be- 



low). Extending the time-dependent solution (28) to the general relativistic case by the prin- 



ciple of minimal coupling (or the convention "comma becomes a semi-colon"), the Michel's 



solution (28) remains valid for arbitrary fi(rfast — t) in General relativity. The argument of 
should be evaluated at the position of a radially propagating fast mode in the Schwarzschild 
metric with driest /dt = o? ^ 

= (r - t + r(l - a^) ^(ra^)) (29) 



The Michel solution in GR has the same flux function as in the fiat space (see Eq. (28)) 



the poloidal magnetic field is derived from $ using a covariant derivative, while toroidal 
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FIG. 1. Cartoon of the structure of magnetic fields around a collapsing rotating neutron star. 
Initially, left panel, the magnetic field is that of an isolated pulsar, with a set of field lines closing 
within the light cylinder (dashed vertical lines). Immediately after the collapse, central panel, the 
structure is similar. The closed field lines are absorbed by the black hole, while the open field lines 
remain attached to the black hole; the system relaxes to the monopole structure (right panel). 

magnetic field and poloidal electric field change according to ^ B^j a and Eq — )■ Ee/a. 
Thus, the exact non-linear general relativistic time-dependent force-free fields corresponding 
to the arbitrary solid-body rotation are 



Br 



Ba 



r 



B 



Si 



-Ds, -c/e 



B, 



-2 



ar 

Is 



r 



a 



(30) 



with n given by Eq. (29). It may be verified by direct calculations that fields (30) satisfy 



the Maxwell equations (12) with the Ohm's law (13) 



As the surface of the neutron star approaches the black hole horizon, Rg — )■ Rg, Bg — )■ 
{Rns/RgYB-^s, while its angular velocity approaches a finite limit which we estimate next. 
Let initially the neutron star rotate with angular velocity ^ns- The moment of inertia of a 
neutron star can be written 

/ns = (2/5)xMNsi?NS (31) 

where x ~ 0.1 — 0.5 is an equation of state dependent variable that describes how centrally 
condensed the star is [35] . The spin angular momentum is thus 



S = (2/5)xMNsi?^s^ 



NS 



(32) 



where Pns is the initial spin period. The dimensionless Kerr parameter is then 



(2/5)x(c/G) 



-^NS^NS 

Mns 
15 



0.04x-iPns,-3 



(33) 



where Pns -3 = -PNs/lmsec. For merging neutron stars the Kerr parameter is expected to 
be much higher. 

For a collapsing star, the time dilation near the horizon and the frame-dragging of the 
horizon lead to the "horizon locking" condition: objects are dragged into corotation with 
the hole's event horizon, which has a frequency associated with it of 

= «2l!^ " f ((S^"''^ = ' 10'-d^"x-.^Nt-3 (34) 

where r^y = (1 + a/1 — a'^)GM/c^ ^ Rq is the coordinate radius of the horizon of the Kerr 
black hole. (Note that for the chosen parameters the final spin is smaller than the initial 
spin, Qh/^tsss = 0.46x-i, due to the assumption of highly centrally concentrated initial 
mass distribution, x ^ 1-) 



The electromagnetic power produced by the Michel's rotator is then (see Eq. (10)) 



2 {BsRI)'^VL'\j 2 gC^-^NS-^NS^H 



It will lead to the black hole spin-down on a time scale 

r = 6 3^',^^^, =2x10^ sec 5^112 (36) 

(Michel solution corresponds to the spin-down index of n = 1, so that the spin evolution 
is described by a decaying exponential.) It is unlikely, though, that the assumptions of the 
model will be applicable for such a long time, see below. 

In addition, the neutron star with dipolar magnetic field has a net charge Q = (1/3)Sns^ns^ 
As long as the assumption of the model are satisfied (that the black hole produces a wind, 
see below), this charge is not canceled by the electrostatic attraction of charges from the 
surrounding medium. Thus, the black hole settles to the Kerr-Newman solution. The 
corresponding Newman parameter is small 



Q^G/c^ Q i?Nsi?^s^NS ^ 4 ^ ^Q_S5^^_^^ p.^_^ (37) 



Rh 2VGM^s 6cMns 

As we argued above, the closed magnetic field lines will be quickly absorbed by the black 
hole, so that the magnetosphere will settle to the monopolar magnetic field structure with 
no electric charge. Fig. [T] 
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B. How a neutron star collapse proceeds 



To summarize the above discussion, first, the space-time of the collapsing neutron star 



temporarily passes through the Kerr-Newman solution with parameters given by (33), (37); 
quickly the electric charge is lost due to the absorption of the closed field lines. (We stress 
that the loss of the electric charge is driven by the internal electrodynamics and not by 
the attraction of charges from the surrounding medium.) Second, and most importantly, 
we have demonstrated that collapsing neutron star does not produce any narrow current 
structures or other dissipative/resistive structure that could have became dissipative and 
"released" the overlaying magnetic field to the infinity: the field always remain connected 
to the surface of the star. 

The fate of the magnetic field lines connected to the surface of the star then depends on 
whether it is a closed magnetic field line, or the one open to infinity. For closed loops, both 
footpoints are dragged toward the horizon and eventually absorbed by the black hole. On 
the other hand, the open magnetic field lines remain open and connected to the hole, without 
"sliding off the black hole", as long as the assumptions of the model remain satisfied. Thus, 
for black hole surrounded by highly conducting plasma the open magnetic field lines never 
become disconnected from the black hole. As a result, the electromagnetic power emitted by 
the black hole may continue for times much longer than the immediate collapse time. 

The key difference here from the conventional BZ mechanism is that in the latter case 
the magnetic field is assumed to be produced by the currents in externally supplied accretion 
disk, while here the magnetic field is produced by the currents generated by the black hole 
itself. Also note, that this result does not violate the "no hair" theorem [e.g., ^36], which 
assumes that the outside is vacuum. In our case the outside medium is assumed to be high 
conducting plasma all the way down to the black hole horizon. Under this assumption the 
magnetic field lines cannot disconnect from the black hole. 

There is a natural limit of applicability of the present model. The electric currents that 
support the magnetic field on the black hole are assumed to be self-produced by the black 
hole via the vacuum breakdown, and not supplied by the external current, like in the BZ 
case. Vacuum breakdown requires a sufficiently high electric potential. As the black hole 
spins down, the potential available for particle acceleration decreases. After some time, 
the black hole will not be able to break vacuum. It would cross a death-line (using pulsar 
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terminology) after which moment no particles are produced anymore, the outside becomes 
vacuum, and by the no hair theorem the black hole will lose its black hole. Also, starting 
this moment the black hole will be able to attract charges of the opposite sign, canceling 
the internal charge. 

In fact, a somewhat different scenario is likely to play out. Our experience with pulsars 
indicate that the plasma production in the magnetosphere is a highly non-stationary process. 
If there is an interruption in the plasma production for sufficiently long time, the magnetic 
field will able to slide off the black hole, shutting down the electromagnetic power forever. 

VII. ON THE NATURE OF SHORT GAMMA RAY BURSTS 

The above results further highlight possible difficulties with the progenitors of short GRBs 
being the merging neutron stars [37] . On the one hand, numerical simulations indicate that 
the active stage of NS-NS coalescence typically takes 10-100 msec. Only small amount, 
< O.IMq of material may be ejected during the merger and accretes on time-scales of 1-10 
sees, depending on the assumed a parameter of the resulting disk [e.g., EBHIO]- Thus, there 
is not enough baryonic matter left outside the BH to power a short GRB. Any energetically 
dominant activity on much longer time scales contradicts the NS-NS coalescence paradigm 
for short GRBs. This seem to contradict observations that some short GRBs have long 
extended X-ray tails observed over time scales of tens to hundreds of seconds. The tail 
fiuence can dominates over the primary burst [by a factor of 30 as in GRB080503, HI]. In 
addition, powerful flares appear late in the afterglows of both short and long GRBs {e.g., 
in case of GRB050724 there is a powerful flare at 10^ sec). In the standard forward shock 
model of afterglows this requires that at the end of the activity, lasting 10-100 msec, the 
source releases more energy than during the prompt emission in a form of low F shells, which 
collides with the forward shock after ~ 10^ dynamical times, a highly fine-tuned scenario. 

On the other hand, the expected electromagnetic powers estimated in the present paper 
are fairly low for all the discussed processes. Since the above results are based on the 
analytical Michel-type structure of the black hole magnetospheres, which for a given surface 
magnetic field and the spin has the largest amount of open magnetic field lines and the 
largest electromagnetic power, the numerical estimates above can be considered as upper 
limits. 
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The only exception to the above could be that an efficient magnetic dynamo mechanism 
operates either during neutron star merger (for short GRBs) or during a core collapse of a 
massive star (for long GRBs), resulting in a formation of a millisecond magnetar-type object 
with magnetic field reaching 10^^ G |12]. Since, as we argue, the black hole can retain its 



magnetic field for a long period of time, the spindown time scale (36 ) may become sufficiently 
short, hundreds to thousands of seconds, so that the magnetic field can electromagnetically 
extract a large fraction of the total rotation energy of the black hole 

^tot = ^Ins^Is = 2 X lO^^ergX-i Pns'-s- (38) 

The fact that the electromagnetic extraction of the rotational energy of the black hole 
can operate both in long and shot GRBs may explain a surprising observation that early 
afterglows of long and short GRBs look surprisingly similar, forming a continuous sequence, 
e.g., , in relative intensity of X-ray afterglows as a function of prompt energy |13]. This 
is surprising in a forward shock model: the properties of the forward shock do depend 
on the external density, while the prompt emission is independent of it. The difference 
between circumburst media densities in Longs (happening in star forming regions) and 
Short (happening in low density galactic or even extragalactic medium) is many orders 
of magnitude. In defense of the forward shock model, one might argue that afterglow 
dynamics depends on Egjecta/n, both of which are orders of magnitude smaller for short 
GRBs if compared with long GRBs. Yet afterglows are very similar and, most importantly, 
form a continuous sequence 

We suggest that the similarity of the early afterglows properties of long and short GRBs, 
,at times < 10^ sec, can be related to the fact that in both cases a spinning black hole can 
retains magnetic field for sufficiently long time to power the prompt and early afterglow 
emission via internal dissipation in the wind 



VIII. DISCUSION 

In this paper we discuss possible electromagnetic signatures of the merging and collapsing 
compact objects. At the in-spiraling stage, in case of NS-NS system, the peak Poynting 
power is Lu^rnax = 6 X lO^^ergs"^ (Sns/IO^^G)^, while for BH-NS systems it is an order of 
magnitude smaller. Both the peak power and the total energy of the precursor emission are 
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fairly small, see [|TTj Only for magnetar-type magnetic fields the corresponding emission can 
be observed at cosmological distances, see [TT] 

We found Michel-type solution for the structure of time-dependent force-free magneto- 
spheres in General relativity. Based on this solution, we argue that contrary to the previous 
estimates the direct emission of the electromagnetic field, powered by the magnetic energy 
stored outside of the neutron star, does not produce a considerable electromagnetic signal: 
such process is prohibited by the high conductivity of the surrounding plasma. 

Most importantly, as long as the black hole is able to produce a highly conducting plasma 
via the vacuum breakdown, magnetic field cannot "slide off" the black hole. As a result, a 
black hole can retain magnetic field for much longer time that is predicted by the "no hair" 
theorem, producing an electromagnetic power for a long time after the collapse, without a 
need for an externally supplied magnetic field. The "no hair" theorem does not apply here 
due to the assumed high conductivity of the plasma surrounding the black hole. (Pulsars 
produce plasma and currents all by themselves, without an external accretion disk.) Since 
in the force-free limit the structures in the current sheet are flying away with the speed 
of light [cf. , the corrugated current current sheet solution in Ref. HI], any magnetic field 
reconnection occurring beyond the light cylinder does not affect the global solution. The 
moment the black hole fails to produce the plasma {e.g., due to spontaneous reconnection 
within the light cylinder), it will quickly lose its magnetic field and stop producing any 
electromagnetic power. (It takes one malfunction to break the black hole electromagnetic 
engine). It will likely to be a random processes, with no typical time-scale, that will terminate 
the EM emission well before the BH spins down. 

I would like to thank Thomas Baumgarte, Vasily Beskin, Scott Hughes, Luis Lehner, 
Harald Pfeiffer, Eric Poisson, Stuart Shapiro and Alexander Tchekhovskoy for discussions. 



[1] V. Ginzburg and L. Ozernoi, JETP 20, 689 (1965). 

[2] T. W. Baumgarte and S. L. Shapiro, ApJ 585, 930 (2003), arXiv:astro-ph/0211339. 

[3] D. M. Eardley and W. H. Press, ARA&A 13, 381 (1975). 

[4] R. N. Henriksen, W. Y. Chau, and K. L. Chan, ApJ 227, 1013 (1979). 

[5] P. D. Morley and I. Schmidt, A&A 384, 899 (2002), arXiv:astro-ph/0201512. 

20 



[6] S. L. Liebling, L. Lehner, D. Neilsen, and C. Palenzuela, Phys. Rev. D 81, 124023 (2010), 
1001.0575. 



[7] 

[8] 
[9 

[lo: 
[11 

[12 

[13 
[14 
[15 
[16 
[17 
[18 

[19 

[20 

[21 
[22; 
[23; 



[24 
[25; 
[26 
[27 
[28; 

[29 

[3o; 



D. A. Uzdensky, ApJ 598, 446 (2003), arXiv:astro-ph/0305288. 
M. Lyutikov, Phys. Rev. D 83, 064001 (2011), 1101.0639. 

S. D. Drell, H. M. Foley, and M. A. Ruderman, J. Geophys. Res. 70, 3131 (1965). 
P. Goldreich and D. Lynden-Bell, ApJ 156, 59 (1969). 

B. M. S. Hansen and M. Lyutikov, MNRAS 322, 695 (2001), arXiv:astro-ph/0003218. 

C. Palenzuela, T. Garrett, L. Lehner, and S. L. Liebhng, Phys. Rev. D 82, 044045 (2010), 
1007.1198. 

S. T. McWilhams and J. Levin, ArXiv e-prints (2011), 1101.1969. 
P. Goldreich and W. H. Juhan, ApJ 157, 869 (1969). 

P. J. Armitage and P. Natarajan, ApJ 567, L9 (2002), arXiv:astro-ph/0201318. 
J. H. Krolik, ApJ 709, 774 (2010), 0911.5711. 

M. Milosavljevic and E. S. Phinncy, ApJ 622, L93 (2005), arXiv:astro-ph/0410343. 

J. R. van Meter, J. H. Wise, M. C. Miller, C. S. Reynolds, J. Centrella, J. G. Baker, W. D. 

Boggs, B. J. Kelly, and S. T. McWilhams, ApJ 711, L89 (2010), 0908.0023. 

P. Chang, L. E. Strubbe, K. Menou, and E. Quataert, MNRAS 407, 2007 (2010), 0906.0825. 

N. Stone and A. Locb, MNRAS 412, 75 (2011), 1004.4833. 

C. Thompson, M. Lyutikov, and S. R. Kulkarni, ApJ 574, 332 (2002), astro-ph/01 10677. 
R. H. Price, Phys. Rev. D 5, 2439 (1972). 

R. D. Blandford, in Lighthouses of the Universe: The Most Luminous Celestial Objects and 
Their Use for Cosmology, edited by M. Gilfanov, R. Sunyeav, & E. Churazov (2002), pp. 

381-+. 

M. Lyutikov, New Journal of Physics 8, 119 (2006), arXiv:astro-ph/0512342. 
A. Gruzinov, ArXiv Astrophysics e-prints (1999), astro-ph/9902288. 
T. Uchida, Phys. Rev. E 56, 2181 (1997). 

S. S. Komissarov, MNRAS 336, 759 (2002), arXiv:astro-ph/0202447. 

K. S. Thorne, R. H. Price, and D. A. MacDonald, Black holes: The membrane paradigm 

(Black Holes: The Membrane Paradigm, 1986). 

X. Zhang, Phys. Rev. D 39, 2933 (1989). 

S. J. Park and E. T. Vishniac, ApJ 337, 78 (1989). 

21 



[31] S. J. Park and E. T. Vishniac, ApJ 347, 684 (1989). 
[32] F. C. Michel, ApJ 180, 207 (1973). 

[33] A. Spitkovsky, ApJ 648, L51 (2006), arXiv:astro-ph/0603147. 

[34] R. D. Blandford and R. L. Znajek, MNRAS 179, 433 (1977). 

[35] G. B. Cook, S. L. Shapiro, and S. A. Teukolsky, ApJ 424, 823 (1994). 

[36] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (San Francisco: W.H. Freeman 
and Co., 1973, 1973). 

[37] M. Lyutikov, in The shocking Universe, edited by G. Chincarini, P. D'Avanzo, R. Margutti, 

and R. Salvaterra (2010), pp. 3-33, arXiv:astro-ph/0911.0349. 
[38] K. Kiuchi, Y. Sekiguchi, M. Shibata, and K. Taniguchi, ArXiv e-prints (2009), 0904.4551. 
[39] Y. T. Liu, S. L. Shapiro, Z. B. Etienne, and K. Taniguchi, Phys. Rev. D 78, 024012 (2008), 

0803.4193. 

[40] J. Faber, Classical and Quantum Gravity 26, 114004 (2009). 

[41] D. A. Parley, B. D. Metzger, J. Granot, N. R. Butler, T. Sakamoto, E. Ramirez-Ruiz, A. J. 

Levan, J. S. Bloom, A. A. Miller, A. Bunker, et al., ApJ 696, 1871 (2009), 0811.1044. 
[42] V. V. Usov, Nature 357, 472 (1992). 

[43] M. Nysewander, A. S. Fruchter, and A. Pe'er, ApJ 701, 824 (2009), 0806.3607. 

[44] S. V. Bogovalov, A&A 349, 1017 (1999), arXiv:astro-ph/9907051. 



Appendix A: Hyperbolic Grad-Shafranov equation 



Generally, we can write 



= -r sin 9V^ xe^ + r sin OQVP xe^ + VF{P) 

8 P 

E = -nVP - —es + BpF' 

w 

E.B = O^F'='^^ (Al) 

The component of the induction equation gives (the poloidal components are satisfied 
identically) 

dtl = ^ (A*F + ro(VQ X VP)) = ^ (A*P + w\B ■ VQ)) (A2) 
The (j) component of the Ampere's law gives the hyperbolic wave Grad-Shafranov equation 
-^^'^' ypj _ Q2p ^ |_4(A*P)/2 - 2IF"dtP{VPf + IF'dt{VPf 
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4(VP • V/)/ + 4:A*Pw^PQ,'^ - 2ro7Q^(VP x V{dtP/Q)) ■ - 2{A*P)IF'dtP 
+ ((VP)2 + 8P)zu^Q{VP ■ Vn) + 8wPn''{WP ■ V(rsin^^))} (yp^2 + 4/2 = ^ 
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